Resolution of sub-wavelength transmission devices formed by a wire medium 
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in [Phys. Rev. B, 73, 033108 (2006)] are studied in this paper using both analytical and numerical 
modeling. It is demonstrated that such transmission devices have sub-wavelength resolution that 
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I. INTRODUCTION 

The resolution of common imaging systems is re- 
stricted by the so-called diffraction limit, since these 
systems operate only with propagating spatial harmon- 
ics emitted by the source. The conventional lenses can 
not transport evanescent harmonics which carry sub- 
wavelength information, since these waves exhibit expo- 
nential decay in natural materials and even in free space. 
Recently, a new type of devices capable of transmitting 
images with sub-wavelength resolution was suggested in 
. These structures are formed by planar slabs of mate- 
rials with specific electromagnetic properties which make 
possible to transfer images from one interface to another 
with sub-wavelength resolution. Such regime of oper- 
ation is called canalization. This regime dramatically 
differs from the sub- wavelength imaging principle formu- 
lated in 0] which employs effects of negative refraction 
and amplification of evanescent waves. The idea is that 
both propagating and evanescent harmonics of a source 
can be transformed into the propagating waves inside of 
a slab of an electromagnetic crystal. Then, these propa- 
gating modes can transmit sub-wavelength images from 
one interface of the slab to the other one. 

In order to implement the canalization regime it is 
required to use materials in which the electromagnetic 
waves have a flat isofrequency contour. Such media are 
capable of transmitting energy only in one direction, al- 
ways with the same phase velocity. The typical example 
of such materials is the wire medium (see Fig. which 
consists of an array of parallel ideally conducting wires 
HLQLUllii; operating at frequencies smaller than its char- 
acteristic plasma frequency. The wire medium supports 
transmission line modes which travel along the wires at 
the speed of light. A slab of wire medium with thickness 
equal to an integer number of half-wavelengths experi- 




FIG. 1: The geometry of the wire medium: a square lattice 
of parallel ideally conducting thin wires. 



ences Fabry-Perot resonance for any angle of incidence, 
including complex ones, and thus such a slab is capable of 
transporting images with sub- wavelength resolution Q. 
The transmission devices formed by wire medium were 
studied both numerically and experimentally in 0. As 
result, the imaging with A/15 resolution has been success- 
fully demonstrated. In |j| it was assumed that the struc- 
ture operates at a very low frequency as compared to its 
plasma frequency and no estimation has been made on 
how the plasma-like properties of the wire medium could 
affect the operation of such devices. In order to reveal the 
restrictions on the resolution caused by plasma-like prop- 
erties of the wire medium, in the present paper we theo- 
retically study the transmission through a wire medium 
slab using the additional boundary condition recently for- 
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mulated in |8(. The wire medium is modeled as a homo- 
geneous material described by an effective permittivity 
tensor. The theoretical results are validated using the 
full wave periodic method of moments |9( , which models 
every detail of the structure and the actual granularity 
of the material. It is proved that the effective material 
model is sufficiently accurate to model the electromag- 
netic response of the wire medium. 

The paper is organized as follows. After a brief in- 
troduction, the expressions of the reflection and trans- 
mission coefficients for a planar slab of wire medium are 
presented. In the third section, the transmission and re- 
flection properties are studied for various frequencies of 
operation using both theoretical and numerical modeling. 
The fourth section is devoted to the study of resolution 
using the Rayleigh criterion. The imaging accuracy and 
the theoretical limit of resolution are studied in the fifth 
section using a half-intensity criterion. Finally, the con- 
clusion is presented. Also, the paper is accompanied by 
an appendix where guided modes in slabs of wire medium 
are investigated. 



II. REFLECTION AND TRANSMISSION 
COEFFICIENTS FOR A WIRE MEDIUM SLAB 

The wire medium is a material characterized by strong 
spatial dispersion even at low frequencies [y]. It can be 
described by a spatially dispersive permittivity tensor of 
the form: 



e = exx + yy + zz, e(u, q x ) = 1 - 



k 2 



medium is also free space. In this paper, we will use this 
result to study the resolution of the transmission devices 
formed by the wire medium. Firstly, we characterize the 
scattering of plane waves by a wire medium slab. 

Let us consider a slab of wire medium with thickness 
d (see Fig. We suppose that the wires stand in free- 
space and are normal to the interface. The wires are 
assumed perfectly conducting, i.e. losses are assumed 
negligible. The distribution of the electromagnetic field 
at the input plane of the transmission device can be de- 
composed in terms of spatial harmonics. The spatial 
harmonics are either propagating plane waves or evanes- 
cent plane waves, and their polarization can be classified 
as transverse electric (TE) or transverse magnetic (TM) 
with respect to the orientation of the wires. In |7| it 
was proved that a slab of the wire medium allows sub- 
wavelength imaging of the part of the spatial spectrum 
with TM-polarization. 
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where the a:- axis is oriented along wires, q x is the x- 
component of wave vector q, k — lo/c is the wave number 
of the host medium, c is the speed of light, and k p = lo p /c 
is the wave number corresponding to the plasma fre- 
quency Up. The plasma frequency depends on the lattice 
period a and on the radius of wires r |lOj| : 
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Note that in this paper we assume that the metallic wires 
are arranged in a square lattice. The permittivity tensor 
is normalized to the permittivity of the host medium. 

The solution of any boundary problem involving the 
wire medium (using an effective medium theory) requires 
the knowledge of an additional boundary condition at 
the interface. In fact, it has long been known that the 
usual boundary conditions (continuity of the tangential 
components of the electroma gnet ic field) are insufficient 
in case of spatial dispersion [Til Il2l Such an ad- 

ditional boundary condition has been formulated in Q 
for the wire medium case. It was proved that the nor- 
mal component of the electric field must be continuous 
at the interface between the wire medium and free space, 
under the condition that the host medium of the wire 



FIG. 2: The wire medium slab. The structure is unbounded 
and periodic along the y- and z-directions. 

To determine the resolution of the imaging, we consider 
that a plane wave with TM polarization impinges on the 
slab. Let H mc be the amplitude of the incident magnetic 
field along the z-direction, and let k = (k x ,k yi Of be 
the wave vector of the incident wave. The x-component 
of the wave vector k x can be expressed in terms of the 
wave number k and of the y-component of the wave vec- 
tor k y as k x = ~j\Jky — k 2 . The incident plane wave 

excites two types of waves in the wire medium: the ex- 
traordinary wave (TM mode) and the transmission line 
(TEM) mode [f| ■ The ordinary wave is not excited since 
it has TE polarization. The tangential component of the 
wave vector k y is preserved at the interface between free 
space and the wire medium. Using this property we can 
evaluate the wave vector associated with each of the ex- 
cited modes in the slab. The wave vector associated with 
transmission line modes is of the form: 



„TEM 



(±k,ky,0) q 



(3) 



These waves travel along the wires with the speed of light 
in free space, independently of the value of transverse 
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component of wave vector k y . On the other hand, the 
extraordinary modes obey to the dispersion equation de- 
rived in 0: 

ql + q 2 y + q 2 z = k 2 - k 2 p , (4) 

and thus the wave vector of an extraordinary mode is of 
the form: 

q™ = (±q x ,k y ,0), q x = -jyfk* + k*-k*. (5) 

If the frequency of operation is smaller then the plasma 
frequency (fc < k p ), then the extraordinary modes are 
evanescent: the x-component of wave vector q x is a purely 
imaginary number for any real k y . In order to emphasize 
this fact it is convenient to use the notation q x = — J7tm 
with 7tm = \J k 2 + ky — k 2 . Also, by analogy, we de- 
note k x — —j'fx with j x = \jk 2 — k 2 . The value of ^ x 

is purely imaginary if k y < k (propagating modes in free 
space), but it becomes a real number if k y > k (evanes- 
cent modes in free space, i.e. the sub-wavelength spatial 
spectrum) . 



The total magnetic field (directed along the z-axis) can 
be written in the following form: 

' e~ jk ' x + Re jk * x , x I) 

£M e -7 TM (i-d/2) 
H(x) _|_^4TM e + 7 TM(a:-d/2) 

~H- = 1 + ^TEM e -ife( K -d/2) 0<x<d (6) 

+J ^TEM e +jk(x-d/2) ^ 

^Te- jk *( x - d \ x>d, 

where R and T are the reflection and transmission coef- 
ficients, A™ and X£ M are the amplitudes the extraor- 
dinary (TM) and transmission-line (TEM) modes in the 
slab, respectively. 

Following H, all components of both magnetic and 
electric fields are continuous at the interface between 
free space and the wire medium. Thus, the magnetic 
field H(x) and both its first dH{x)/dx and second 
d 2 H(x)/dx 2 derivatives by x are continuous at interfaces 
x = and x — d. This can be written using H in the 
form of the following system of equations: 
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Solving this system the reflection and transmission coefficients can be expressed as follows: 
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T = - - - (9) 

1+ Tl"fL tanh(^) J^-^tan^) 1+ TI^L ctanh(^) + ^ ctanfg) 

7x(feJ+fe2) V. 2 > 7^(fey+fc|) K 2 ' 7x(fey+fep) \ 2 > ' 7x(fcJ+fc^) v 2 > 



In the above "tanh" and "ctanh" represent the hyper- 
bolic tangent and cotangent, respectively. The canaliza- 
tion regime is observed for the slab under consideration 
when kd = nir, where n is an integer. In this case T can 
be simplified: if kd = (2n + l)ir then 



T = — 



1 + 7x ]™+^) ctanh( 7T M^/2) 



R = 1 + T, (10) 



whereas if kd 
T = 



2nir then 
1 



1 



7TMfcg 
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tanh(7TM^/2) 



R=l-T. (11) 



The transmission coefficient calculated above can be 
regarded as a transfer function: the image at the output 
plane is obtained by superimposing the (TM-polarized) 
spatial harmonics of the source field distribution at the 
input plane multiplied by T. Strictly speaking, that is 
not rigourously true, even if the exact expression for 
transmission coefficient is known. Indeed, when doing 
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this we are neglecting the higher order harmonics, which 
decay much faster than the fundamental mode, but also 
contribute to the near field. These harmonics are a man- 
ifestation of the actual granularity of wire medium, and 
are not described by the effective medium model eff. The 
contribution of the high order harmonics can be neglected 
while k y a < tt, which corresponds to operation below the 
ultimate limit of resolution A; = a 14] . In this case, 
slight variations of the near field near the edges of the 
wires do not influence the resolution of the device. Thus, 
in this paper we will consider that the distribution of the 
field at the output plane is equal to the distribution of 
the field at the input plane multiplied by the transmission 
coefficient. Within this approximation, it is sufficient to 
characterize T in order to evaluate the resolution of the 
transmission device. 

If the operating frequency is much smaller than the 
plasma frequency, i.e. k <C k p , then since kd/n > 1 we 
can assume that 7tm^ 3> 1, and using the expressions 
t2nl and t2n we obtain simple approximate formulae for 
the reflection and transmission coefficients: 

T~T 1 , 2 , iZwl 1 , (12) 

where the =p sign corresponds to the case in which the 
thickness of the slab is equal to an even or odd number 
of half- wavelengths, respectively. 

Some important properties of t2nl and t2n are readily 
identified. It is clear that if k p — > oo then T w =pl for any 
k y , and therefore, in such conditions, the imaging is per- 
fect for this polarization. On the other hand, one can see 
that T — when j x = 0. i.e. k y = k. This property is 
very important because it shows that the spatial harmon- 
ics with k y = k are filtered by the WM, which some how 
contradicts the sub-wavelength imaging. Actually, as it 
will be demonstrated below, the band of the spatial har- 
monics corresponding to this filtering is so narrow that it 
practically does not affect the imaging properties of the 
system. 

It is noteworthy that if k y — > oo then T — > =Fl/2 [2lj| . 
This means that the spatial harmonics in the deep sub- 
wavelength spectrum are transferred from the source to 
the image plane, but in contrast to the spatial harmon- 
ics with smaller k y , their amplitude will be reduced two 
times. Thus, the transfer function of the imaging device 
is not perfect since it is not a constant function of the 
transverse wave vector. Nevertheless, it will be shown 
ahead that this fact does not significantly affect the sub- 
wavelength imaging performance of the device. 

Also, it is clear that the imaging quality does not de- 
pend on the thickness of the structure. The transmission 
device can be made as thick as it may be required by 
an application. The only restriction is that the thickness 
has to be equal to an integer number of half- wavelengths. 
If the number of half-wavelengths is odd then the image 
at the back interface of the structure will appear out-of 
phase with respect to the source, whereas if the number 
of half- wavelengths is even then the image and the source 



will be in phase. Moreover, the bandwidth of operation 
does not depend on the thickness of the transmission de- 
vice. It remains the same for all slab thicknesses equal 
to the integer number of half- wavelengths. 



III. PARAMETRIC STUDY OF REFLECTION 
AND TRANSMISSION PROPERTIES 

In the wire medium with period a — 1 cm formed 
by wires with r = 1 mm has been considered. Substitu- 
tion of these parameters into kO gives k p a = 2.5. A more 
accurate value, calculated using the numerical method 
proposed in [15|. is k p a = 2.36. The frequency of opera- 
tion for the transmission device in Q was 1GHz, which 
corresponds to ka = 0.2, and the thickness d = 15 cm of 
the structure was chosen equal to half-wavelength, that 
is kd = 7r, or in terms of the wave number corresponding 
to the plasma frequency k p d/ir — 11.3. In the present pa- 
per we consider the same parameters for the transmission 
device. 
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Transverse component of wave vector, k y d/n 

FIG. 3: (Color online) Reflection (R) and transmission (T) 
coefficients (absolute values) as a function of the transverse 
component of wave vector k y d/n for kd/n = 1. Solid lines: 
analytical expressions R and T. Crosses and circles: points 
calculated using full wave numerical method. 

The dependence of the reflection and transmission co- 
efficients (absolute values) on the normalized transverse 
component of wave vector k y d/ir is plotted in Fig. [3] 
One can see that the transmission coefficient has a very 
sharp dip at k y = k, but for other values of k y it has 
values close to unity. In Fig. [3] the "crosses" and "cir- 
cles" represent the data calculated numerically using the 
periodic method of moments. As seen, the agreement 
between the analytical model and the full wave results is 
very satisfactory. We do not present a plot for the phase 
of transmission coefficient since it is practically equal to 
7r for whole range of k y , except for a very narrow band 
with upper bound equal to k. The obtained dependence 
confirms that the slab indeed operates in the canalization 
regime and is capable of transporting sub- wavelength im- 



ages with TM polarization. 
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Transverse component of wave vector, k y d/n 

FIG. 4: (Color online) Reflection (R) and transmission (T) 
coefficients (absolute values) as a function of the transverse 
component of wave vector k y d/n for kd/n = 1 — 0.002n, where 
n=0,l,2,3,4. The numbers in the figure correspond to the 
values of kd/n. 

The behavior of the transmission and reflection coeffi- 
cients is sensitive to variations in the frequency of opera- 
tion. In Fig. 01 the transmission characteristic is depicted 
for frequencies that are lower than original frequency of 
operation by 0.2,0.4,0.6 and 0.8%. It is seen that the 
amplitude of the evanescent modes is amplified in a cer- 
tain range of k y . 

The further decrease of the frequency of operation re- 
veals a very interesting resonant phenomenon. This is 
illustrated in Fig. \5\ where the transmission coefficient 
is plotted for frequencies which lower than the design 
frequency by 1,2,3,4 and 6%. One can observe a very 
strong enhancement of certain spatial harmonics. This 
property reveals the presence of guided modes propagat- 
ing along the y direction of the slab, being the mechanism 
of pro pag ation closely related to that of an Yagi antenna 
array [16| . The dispersion properties of the guided modes 
are studied in the Appendix. The observed resonant en- 
hancement deteriorates the sub- wavelength imaging since 
some of the spatial harmonics happen to be amplified by 
a factor of 10 or even more. 

For frequencies slightly larger than the frequency cor- 
responding to the Fabry-Perot resonance, the reflection 
and transmission coefficients do not have a resonant be- 
havior. This is shown in Fig. where the reflection and 
transmission coefficients arc plotted for frequencies larger 
than the design frequency by 1, 2 and 4%. This means 
that the slab does not support guided modes in this fre- 
quency range. Further increase of the frequency (by 8 or 
16%, for example, see Fig. EJ) leads to the increase of the 
width of the narrow dip of the transmission character- 
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FIG. 5: (Color online) Reflection (R) and transmission (T) 
coefficients (absolute values) as a function of the transverse 
component of wave vector k y d/n for kd/n = 1 — O.Oln, where 
n=0,l,2,3,4,6. The numbers in the figure correspond to the 
values of kd/n. The crosses and circles represent the results 
of the numerical modeling for kd/n — 0.94 and kd/n = 0.96. 
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Transverse component of wave vector, k y d/n 

FIG. 6: (Color online) Reflection (R) and transmission (T) 
coefficients (absolute values) as a function of the transverse 
component of wave vector k y d/n for kd/n — 1 + O.Oln, where 
n=0,l,2,4,8,16. The numbers in the figure correspond to the 
values of kd/n. The crosses and circles present the results of 
the numerical modeling for kd/n = 1.04 and kd/n — 1.08. 

istic at k y — fc, and to the appearance of resonances in 
the vicinity of this point. These resonances (shown with 
greater detail in Fig. [7J are also related with guided 
modes. The dispersion branch of these modes appears 
near the light line in Fig. The resonant behavior in 
the present case has a very narrow band character and 
does not significantly affect the imaging properties of the 
device. 

In Fig. one can see a strong reduction of the reflec- 
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Transverse component of wave vector, k y d/n 



FIG. 7: (Color online) Reflection (R) and transmission (T) 
coefficients (absolute values) as a function of the transverse 
component of wave vector k y d/n for kd/n = l + 0.04n, where 
n=0,l,2,3,4. The numbers in the figure correspond to the 
values of kd/n 

tion coefficient for kd/ir = 1.04 and kd/ir — 1.08. The 
reflection coefficient becomes less than 10% practically 
for the whole spatial spectrum. This fact has been nu- 
merically confirmed, and means that for these frequencies 
there is little interference between the source and the sig- 
nal reflected at the wire medium slab. 



IV. STUDY OF RESOLUTION 

In order to study the resolution of the wire medium 
slab, we use the Rayleigh criterion: the resolution is taken 
equal to the radius of the image spot (at the half- intensity 
level) for a very sharp (nearly point) source. The resolu- 
tion value obtained using this criterion corresponds to the 
half of the minimum distance between two sharp maxima 
that can be resolved using the imaging device. 

In order that the wire medium homogenization model 
can be used we are not allowed to consider any sources 
with radius much smaller than the real period of the 
structure. Indeed the spatial spectrum of such a source 
contains important spatial harmonics in the range k y a 
7r, whose propagation can not be described by homoge- 
nized model. Moreover, it was shown in [14| that the 
resolution of any periodic structure is limited by its pe- 
riod. Due to these reasons, below we consider a magnetic 
field distribution with diameter only slightly smaller than 
the period. The spectrum of the source at the front plane 
of the wire medium slab is, 

s(k y ) = e -\A?- fe2a / 2 . (13) 

Note that the spectrum is independent of fc z , and con- 



sequently the field intensity is uniform along the z- 
direction. This two-dimensional source has the magnetic 
field distribution: 

H z (y)=A f s(k y )e-^yydk y (14) 



= AjnH™ Ma/2) 2 + y 2 ) = , 

v ' v( a / 2 ) +y 

(2) 

where H{ (x) is the Hankel function of the second kind, 
and A is a constant that defines the amplitude of the 
source. The intensity of the field normalized by its max- 
imum value H z iax = H z (0) is plotted in Fig. |H1 (thick 
solid line). 

The distribution of magnetic field under consideration 
appears at the half-period distance from a plane with 
point source of magnetic field which has uniform spatial 
spectrum. Thus, expression spectrum also represents the 
transmission coefficient for a slab of free space with the 
half-period thickness. This slab allows to transform a 
singular distribution of the point source into less sharp 
distribution which is appropriate for our studies using 
effective medium theory. Note, that the point source of 
magnetic field we are dealing with has nothing to do with 
magnetic line source (line of magnetic current flowing 
along z-direction)! The point source of magnetic field is 
a delta function of magnetic field in the plane x = —a/2 
[H z (y)\ x= _ a /2 — S(y)] which means that the field in that 
plane is zero everywhere except at one singular point. We 
have chosen this source because of the simplicity of its 
exponentially decaying spatial spectrum spectrum and 
since its field decays faster than the field created by a 
point magnetic source: the field of our source source de- 
cays inversely proportional to y, but the field radiated by 
a magnetic line source would decay as the inverse square 
root of y. 

The magnetic field at the output plane when the slab 
is illuminated by the source is given by: 

H z {y) = A f s{k y )T(k y )e-^y y dk y . (15) 



The integral in image is singular and non-integrable 
if the slab supports guided modes: some poles of the 
transmission coefficient T as a function of k y are located 
at the real line (see the Appendix). These difficulties 
can be avoided if negligibly small losses are introduced 
into the permittivity of the host medium. The presence 
of losses shifts the poles away from the real line. This 
makes the integral non-singular and it can be relatively 
easy evaluated using direct numerical integration. The 
integral image has been calculated for all values of kd/ir 
considered in the previous section. The corresponding 
images at the output plane are plotted in Fig. [H] and 



7 




Transverse coordinate, y/a 



FIG. 8: (Color online) Image when the wire medium slab 
is illuminated with the source source at frequencies slightly 
smaller than kd/ir = 1. The numbers in the figure correspond 
to the values of the normalized frequency kdjix. 
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Transverse coordinate, y/a 

FIG. 9: (Color online) Image when the wire medium slab 
is illuminated with the source source at frequencies slightly 
larger than kd/ir — 1. The numbers in the figure correspond 
to the values of the normalized frequency kd/ir. 

[5] for frequencies below and above the first Fabry-Perot 
resonance, respectively. For ka < 0.93 and ka > 1.5 the 
image is distorted by the resonant guided modes of the 
slab with short and long wavelengths, respectively. On 
the other hand, for 0.93 < ka/ir < 1.5 the images have 
the same shape as the source and their radii vary from a 
to 2a. This means that at these frequencies the resolu- 
tion of the imaging device A is from A/30 to A/15. Note, 
that A; = A/30 = a is the ultimate limit of resolution 
imposed by the periodicity of the structure 14]. The 
resolution of A/15 for the structure under consideration 
has been observed in both numerically and experi- 
mentally. This allows us to conclude that the theoretical 
study presented in this paper is consistent with results 
reported in 0. However, the presented analysis gives a 
much more comprehensive explanation of the limitations 



and characteristics of our imaging device. 



V. ACCURACY OF IMAGING 

The Rayleigh criterion used in the previous section 
is classical for imaging above diffraction limit. How- 
ever, near-field sub-wavelength imaging significantly dif- 
fers from conventional imaging. Indeed, the latter re- 
quires that the source is placed very close to the interface 
of the imaging device. Hence, not only the transmission 
properties of the structure affect its performance, but also 
the reflection properties. If the reflection coefficient has 
large amplitude then the reflected field can interact with 
the source and modify its near field pattern. This modi- 
fied pattern will be transmitted with sub-wavelength res- 
olution from the input plane to the output plane, but 
then it has nothing to do with the original source. As 
discussed in the previous section and seen in Figs. 
and|Hl the strong reflections from the slab of wire medium 
are mainly caused by the resonances associated with the 
guided modes. That is why the imaging with minimum 
reflections happens at the frequencies corresponding to 
the band gaps for the guided modes and slightly higher 
frequencies. Fig. clearly illustrates this statement: for 
1.01 < kd/ir < 1.16 the reflection coefficient from the 
slab is smaller than 20%, and for for 1.04 < kd/ir < 1.08 
it is smaller than 10%. 

Another important point concerning accuracy of the 
sub- wavelength imaging using slabs of wire media is that 
the transmission coefficient is different for different trans- 
verse components of the wave vector. As it was already 
mentioned, even if the thickness of the slab is tuned to 
fulfil Fabry-Perot condition (as in Fig. |3J) the transmis- 
sion coefficient is close to unity only in a limited range 
of spatial harmonics. For spatial harmonics with large 
k y the transmission coefficient approaches 1/2. This fact 
causes distortion in the image. It is possible to estimate 
the resolution at which the imaging device works with 
acceptable distortion. With this purpose we use the a 
half-intensity criterion: we assume that imaging with ac- 
ceptable distortion happens while the transmission coef- 
ficient of spatial harmonics is in the range [l/y/2, y/2]. If 
fc™ ax is the maximum transverse wave vector component 
for which the half-intensity criterion is fulfilled then the 
resolution is A = 7r/fc™ ax . This definition of resolution is 
stronger than the Rayleigh criterion. It guarantees not 
only that maxima at distance A are resolved, but also 
that the image is restored with little distortion. 

For the case presented in Fig. [21 the transmission 
coefficient T is less then 1 and greater than l/v2 for 
kyd/ir < 5 except for the very narrow range of k y d/ir 
close to 1 . This means that the resolution of the transmis- 
sion device under study is equal to A = 7r/fc," lax = A/10. 
This value is only three times larger than the ultimate 
limit of resolution due to periodicity of the structure 
A; = a = A/30 formulated in 0], and a bit worse than 
the real resolution of 2a = A/15 observed at 0. 
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A slight change of frequency may help improving the 
resolution of the device. For example, in the case of 
kd/n — 0.996 (see Fig. 0J) the resolution happens to 
be equal to A/34 which is even smaller than the ulti- 
mate limit dictated by periodicity. When kd/ir = 0.96 
the resolution is reduced down to A/10. One can ver- 
ify that imaging with a resolution better than A/ 10 is 
observed for kd/n G [0.96,1.00]. This means that the 
bandwidth of operation with A/10 or better resolution 
is equal to 4% for the transmission device under con- 
sideration. The bandwidth of imaging with worse but 
still sub-wavelength resolution can be even larger. In Q 
the 15% bandwidth of imaging with sub-wavelength res- 
olution has been reported. We summarize our results 
for the dependence of the resolution (defined using the 
half-intensity criterion) on the frequency in Fig. 1101 By 
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FIG. 10: (Color online) Normalized resolution A/A as a func- 
tion of the normalized frequency of operation kd/n. 

tuning the frequency one can dramatically enhance the 
resolution. In practice this means that it is always pos- 
sible to reach an ultimate limiting resolution A; = a 
|14| by appropriately choosing the frequency of opera- 
tion. The value of the resolution for the case kd/n = 1 
approximately describes an average level of resolution of 
the system. This value can be evaluated analytically as 
explained next. 

From equation app it is clear that for k y k (sub- 
wavelength spatial spectrum) then 



1 



(16) 



Solving the equation T(fc™ ax ) = 1/V2 we obtain 



t max 



'V2-1 



0.455fc B 



(17) 




5 10 15 20 25 30 35 
Transverse component of wave vector, k y d/n 

FIG. 11: (Color online) Dependence of the transmission co- 
efficient T (absolute value) on transverse component of wave 
vector kyd/n for the case when kd/n = 1 and various values of 
plasma frequency k p d/n — 12.5 ■ 2 n , where n=0,l,2,3,4. The 
numbers at the figure correspond to the values of k p d/n 



by: 



Thus, the resolution A = 7r/fc™ ax for kd/n = 1 is given 



A = 



0.455fc o 



= 1.1A— = 2.75a* /In 



2nr 



0.5275. 



(18) 

Thus, the smaller the ratio k/k p is, the better is the res- 
olution of the system. Formula fesol shows that the res- 
olution of the wire medium slab is only limited by the 
ability to fabricate very dense arrays of wires, i.e. the 
limit of resolution only depends ultimately on the value 
of the lattice constant of the crystal. By decreasing the 
period of wire medium one can greatly improve resolution 
of the system. This fact is illustrated in Fig. ^] which 
shows the transmission characteristic calculated for dif- 
ferent values of the the plasma frequency. For extremely 
thin wires and very high frequencies the effect of losses 
may not be negligible. Nevertheless, at the microwave 
domain, this effect is expected to be of second order. 



VI. CONCLUSION 

In this paper, the resolution of sub-wavelength trans- 
mission devices formed by the slabs of wire media was 
studied. It was shown that the resolution does not de- 
pend on the thickness of the structure, which can in prin- 
ciple be made as thick as required by an application. By 
slightly tuning the frequency of operation, it is possible 
to achieve the ultimate limit on resolution dictated by the 
periodicity of the system (lattice constant of the crystal) , 
even though this regime is narrow band. The average 
level of resolution ultimately depends on the period of 
the lattice. By reducing the lattice constant it may be 



9 



possible to realize imaging systems with virtually no limit 
of resolution. This makes the slab of the wire medium a 
unique imaging device capable of transmitting a distri- 
bution of TM-polarized electric field with sub- wavelength 
resolution at the microwave frequency range. 
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APPENDIX A: GUIDED MODES IN THE SLAB 
OF WIRE MEDIUM 



The dispersion equation disp is transcendental, and so 
cannot be solved analytically. Its numerical solution is 
presented in Fig. 1121 in the form of a dispersion dia- 
gram. The parameters of the wire medium slab are the 
same as in the whole paper: a = 1 cm, r = 1 mm and 
d = 15 cm (kpd/ir = 11.3). One can see that the slab 
of wire medium supports guided modes (when q is a real 
number and q > k) nearly for all frequencies, except for 
very narrow bands that occur near resonances that cor- 
respond to an integer number of half-wavelengths across 
the slab [d = n.A/2, where n is an integer). In this 
case the eigenmodes are leaky waves (Im(g) ^ and 
Rc{g} < k). The dispersion curves for the leaky modes 
are presented in the shadowed region in Fig. Only 
the real part of the propagation constant q is shown for 
the leaky modes. 



A closely-spaced chain of half-wavelength wires 0] 
(one-dimensional Yagi antenna array) can be regarded as 
a waveguide formed by resonant scatterers. This struc- 
ture is in a certain sense the microwave analogue of op- 
tical plasmonic waveguides [ill ITsl IT^ | . A slab of wire 
medium can be considered as a two-dimensional analogue 
of the one-dimensional Yagi antenna array. The eigen- 
modes of such a waveguide were studied in j2lj, but the 
wire medium was modeled as a local dielectric. Below, we 
derive the dispersion equation for the eigenmodes taking 
into account both the spatial dispersion effects Jo], and 
the additional boundary condition proposed in 

Let us consider a slab of wire medium with thickness 
d (see Fig. |2Jl, and investigate if this structure can sup- 
port guided modes travelling along the y-direction with 
some propagation constant q. A guided mode must have 
TM polarization with respect to the wires, since the wire 
medium is transparent to the other polarization. Then, 
the total magnetic field (directed along the z-axis) can 
be written as in H, but without the term corresponding 
to the incident wave since the guided mode is only sup- 
ported by its own field. The application of the boundary 
conditions formulated in y| (the continuity of all com- 
ponents of both magnetic and electric fields at the inter- 
face between free space and the wire medium and, con- 
sequently, the continuity of the magnetic field and of its 
both first and second derivatives by x) at both interfaces 
x = and x = d provides a homogeneous system of equa- 
tions similar to syst (the system is homogeneous due to 
the absence of the incident wave). A non-trivial solution 
of the system exists only if the determinant of the ma- 
trix in syst vanishes. This yields the following dispersion 
equation: 



1 - 



7 TM^tanh(^) - kkl tan (f)' 



1 + 



7TM^ctanh ( 



7tm^ 
2 



+ kkpCt&n ( 



2 / 



7*(*2 + *p) 



= 0. 
(Al) 
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FIG. 12: (Color online) Dispersion diagram for the wire 
medium waveguide. The region above the light line (q — k) 
corresponding to the leaky wave region is shadowed. 



The peaks of the transmission coefficient observed in 
Figs. \5\ and [7| are caused by the resonant excitation of 
the guided modes described in this Appendix. Math- 
ematically, this effect can be readily explained, because 
when the dispersion equation disp is fulfilled the determi- 
nant of the matrix in syst vanishes, and consequently the 
reflection and transmission coefficients R, T have poles. 
Thus, when an evanescent plane wave with transverse 
component of the wave vector k y equal to the propaga- 
tion constant q of the guided mode illuminates the slab, 
the transmitted wave is very much amplified because of 
the described resonant phenomenon. Note, that in Fig. 
the resonances are absent for 1 < kd/w < 1.06 since the 
guided modes do not exist in the considered frequency 
range and the leaky modes cannot be excited. 
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